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Strain-induced gradients of local electric fields in semiconductor quantum dots can couple to 
the quadrupole moments of nuclear spins. We develop a theory describing the influence of this 
quadrupolar coupling (QC) on the spin correlators of electron and hole "central" spins localized in 
such dots. We show that when the QC strength is comparable to or larger than the hyperfine coupling 
strength between nuclei and the central spin, the relaxation rate of the central spin is strongly 
enhanced and can be exponential. We demonstrate a good agreement with recent experiments on 
spin relaxation in hole-doped (In,Ga)As self- assembled quantum dots. 

PACS numbers: 



The spin of an electron or a hole in a semiconductor 
quantum dot is the main component of numerous pro- 
posed spintronic and quantum computing devices^. Spin 
decoherence and finite spin lifetimes are currently the 
major factors that limit our ability to control spin states 
in dots. A single "central" (i.e., electron or hole) spin in 
a dot interacts via hyperfine coupling with a large num- 
ber (10 4 — 10 6 ) of nuclear spins. The net effect of this 
coupling to the nuclear spin bath can be characterized 
by an effective Overhauser magnetic field B n that acts 
upon the central spin. Within a quantum dot ensemble, 
each central spin precesses around a different B n . If B n 
is time-independent, such precession alone cannot lead 
to complete relaxation of the central spin polarization. 
This is evidenced from the observation of spin echoes 2 
that can be used to cancel the dephasing of central spins 
in an ensemble of dots with different constant B n . How- 
ever, stochastic dynamics of the Overhauser field B n in- 
duces irreversible relaxation of the central spin and loss 
of coherence^'. The physics that leads to changes of B n 
and its corresponding influence on central spin relaxation 
are the subject of considerable theoretical debate^H^. 

It was suggested that, at microsecond time scales, 
the dynamics of the Overhauser field is dominated by 
hyperfine-mediated nuclear co-flips, which originate from 
unequal strengths of the hyperfine couplings of the cen- 
tral spin to different nuclear spins inside the same dot 4 . 
Numerical simulations by Al-Hassanieh et alP showed 
that such co-flips generally lead only to a logarithmically 
slow central spin relaxation. In contrast, recent experi- 
mental studies with hole-doped (In,Ga)As quantum dots 
reported a nearly ideal Lorentzian shape of the spin noise 
power spectrum, indicating exponential relaxation of cen- 
tral hole spins rather than a power-law or logarithmic 
relaxation 10 . 

Here we show that quadrupolar couplings (QC) of nu- 
clear spins to the strain induced electric field gradients 
inside typical semiconductor quantum dots can induce 
relatively fast dynamics of the Overhauser field B n , and 
consequently accelerated relaxation of electron and hole 
spins in weak external fields. Our model directly applies 



to InGaAs self-assembled quantum dot systems, which 
are among the most popular platforms for spin memo- 
ries and qubitd^^l however, the model applies gener- 
ally to all dots composed of quadrupolar-active nuclei. 
We model such a nuclear spin bath by introducing static 
fields acting on nuclear spins due to QC, in addition to 
the hyperfine couplings to the central spin. We numer- 
ically compute the dynamics of our model by applying 
a time-dependent mean field (TDMF) algorithm^ that 
allows us to study the relaxation of a central spin cou- 
pled to an unpolarized spin bath containing up to ten 
thousand nuclear spins. 

At low temperatures and at time scales shorter than 
a millisecond, a Hamiltonian that captures central spin 
dynamics in quantum dots has the following form: 

N 

H = Y^ {l\\IizS z + t1(4A + 4A)) + 9zB z S z + 

i=l 

9X BJ X + gyBySy + £ ^ ^(^ • 71^ ~ , (l) 

1=1 ^ ' 

where S and Ii stand for spin operators of, respectively, 
central and nuclear spins; B a is an applied magnetic 
field component along the axis a; g a is the correspond- 
ing component of the central spin g-factor. Index i runs 
though all nuclear spins that interact with the central 
spin. Parameters 7n and 7^ are the out-of-plane (lon- 
gitudinal) and in-plane (transverse) coupling strengths, 
respectively, between the central spin and z-th nuclear 
spin. Henceforth we drop index i for coupling strengths 
when we discuss their typical magnitudes. For electrons, 
71 1 and 7^ have similar magnitudes, but 711 and 7^ are 
quite different for holes. For the latter case, the ratio of 
transverse to out-of-plane couplings, f3 = 7_l/7|| varies in 
different samples in the range 10 13 (3 ~ 0.1 — 0.7. Addi- 
tional coupling terms in the Hamiltonian such as ~ S z Ii X 
are allowed but they were estimated to be negligibly small 
both in electron and in hole-doped dots^l, and we will 
disregard them. We also disregard the Zeeman coupling 
between the external field and nuclear spins because we 
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consider only weak external fields, about the size of the 
Overhauser field (~ 25 Gauss for an InGaAs hole doped 
doPj). 

The last term in ([9| describes QC with characteris- 
tic strength 7q, and coupling anisotropy vector for 
the i-th nuclear spin. QC is allowed for nuclei having 
spin larger than 1/2. QC has previously proved impor- 
tant in experiments on polarized spin bath relaxation in 
GaAs^Mill however, it has been generally disregarded in 
the context of central spin relaxation with initially un- 
polarized nuclear spin baths, both in electrorP^ and in 
hole-doped^ dots. We believe that this omission can- 
not be justified except in certain materials, such as Si, 
that contain predominantly spin-0 or spin- 1/2 nuclei. In 
the widely studied InGaAs dot system, the most abun- 
dant indium isotopes 115 In and 113 In have / = 9/2, and 
Ga an d As iso topes have / = 3/2. According to many 
7q ~ 2 — 4 MHz for indium atoms in 
GaAs at a typical strain of 3-4% inside a dot, which 
translates for spin 9/2 into a characteristic level split- 
ting 7 C = Jq\I\ ~ 10 MHz. This value is at least 
an order of magnitude larger than the effective hyper- 
fine coupling 7|| ~ 0.1-0.5 MHz in a typical hole-doped 
quantum dot with N - 10 5 nuclei 1 ^. Recent NMR 
studies of InGaAs dots also showed that the directions 
of QC anisotropy axes are strongly non-uniform in- 
side a dot and do not align with the sample growth 
anisotrop}ff^. To include this fact, we will assume that 
the local anisotropy vector rii for the i-th nuclear spin 
points in a random direction, which is chosen indepen- 
dently for each nuclear spin. Note that this does not ex- 
clude arbitrary spatial correlations of different inside 
the dot. 

The Hamiltonian (|9| belongs to the class of spin bath 
models, in which noncollinear static fields act on nuclear 
spins independently of the coupling to the central spin. 
In order to compare different models of this class, we 
introduce a parameter j c that characterizes the typical 
energy level splitting of nuclear spins by static fields. Our 
theory shows that this parameter determines all the es- 
sential effects of the static fields irrespective of the details 
of the interactions. This renders our theory applicable to 
spin baths with different sizes of nuclear spins /. Our 
results extend beyond the Hamiltonian In fact, the 
minimal model of our class of spin baths can be formu- 
lated in terms of the central spin problem with only a 
nuclear spin- 1/2 Hamiltonian: 



N 



(2) 



where the last term mimics the effect of QC, & is the 
Pauli operator of the central (electron or hole) spin, & % a 
is the a-component of the Pauli operator for the i-th nu- 
clear spin, and 7* corresponds to the size of the charac- 
teristic level splitting for the i-th nuclear spin with quan- 
tization axis n\ 




FIG. 1: Spin correlator C 2 (t) = (a z (t)a z (0)) at -y c = 0, 7 N = 
(7ll) = 1, N = 700 nuclear spins, and B = 0, shown up to 
times (a) t = and (b) t — 0.5/7|| (time t is in units of 

1/7||). Here /3 = 7±/7|| is the hyperfine coupling anisotropy. 



We will compare, in Fig. |2j the dynamics of the model 
with the Hamiltonian ([9| for the spin bath with 1 = 1 
and the minimal model ([5| at the same characteristic 
value of 7 C . Results are almost indistinguishable, so in 
the rest of the main text, we will show numerical results 
only for the minimal model to illustrate all the effects. 

In the supplementary file 21 , we describe the TDMF ap- 
proach and provide additional numerical tests for evolu- 
tion with N from 250 to 10000 nuclear spins, the Hamil- 
tonian ([9| with 7=1, and the classical limit I ^> 1, 
which are all found to be in very good agreement with 
the theory that we develop here. 

Here, we present our results for central spin tempo- 
ral correlators that were obtained for the model ([5| with 
N = 700 spin- 1/2 nuclei at equilibrium. Before each 
simulation, we chose 7^ = 2j\\ * r^, 7^ = 2/3j\\ * T2% and 

7* = 27 c *r3^, where r^, r^i and r^i are random numbers 
from a uniform distribution in the interval (0, 1). We set 
the energy scale so that 7|| = 1. Note that we chose 
widths of parameter distributions to be comparable to 
the mean values as suggested in pQ. Vectors n % point in 
random directions and the time step was dt = 0.0001. 
Averaging was performed over 1000 and over 30000 ran- 
domly chosen initial state vectors (both for central and 
nuclear spins) for the calculation of, respectively, noise 
power and real-time correlators. 

In Fig. 13 'a) we show our numerical results for the 

(cr z (t)d- z (0)}, obtained 



central spin correlator, C^t) 

from the evolution of the Hamiltonian ([5| in the absence 
of quadrupolar interactions (j c = 0). Different curves 
correspond to different values of the coupling anisotropy 
(3. All curves start at 6^(0) = 1. Figure [l3]fo) resolves 
the part of Fig. [l3^a) with t < 0.5. The appearance and 
shape of the deep local minimum of C2(t) in Fig. 13 'b) is 
well understoocp^ as being due to dephasing caused by 
ensemble central spin precession around the Overhauser 

fields B n = Eti[7i(4)x + 7l(4)y + 7f, <*£>*]• Fi S~ 
ure [l3^a) shows that a fraction of the central spin po- 
larization additionally relaxes during a longer time in- 
terval that is of order 1 / (71 1 ) . This relaxation follows 
from the co-flip effect 1 56 . Figure. 13 'a) confirms previ- 
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FIG. 2: (a) The real time central spin-spin correlator for dif- 
ferent magnitudes of the static field j c in spin- 1/2 bath, (b) 
Exponential fit (dashed red) of the spin correlator (blue) for 
spin- 1/2 bath at j c = 12 in units of 7||. (c) Central spin 
correlator in spin bath with 7=1 and quadrupole coupling 
7q = 4j l c . (d) Exponential fit (dashed red) of spin correlator 
in spin-1 bath at 7 C = 12. In all cases: f3 — 0.2 and B — 0. 



oils observation^, which was made for the case /3 = 1, 
that only a fraction of the central spin polarization re- 
laxes via this mechanism on time scales of interest. It also 
shows that the correlator decay is strongly suppressed by 
hyperfine coupling anisotropy. This means that standard 
co-flip effect cannot explain the observed spin relaxation 
at a fracti on of a m icrosecond in experiments with hole- 
doped dot d 1Q | n |2 ^, for which 1/(7||) corresponds to sev- 
eral microseconds^. 

Figure [2] shows the central spin correlator C2(t) for sev- 
eral different mean values of QC [as tuned by the static 
field 7 C to compare spin-1/2 bath in Figs.[2ja-b) and spin- 
1 bath with Jq = 4j l c in Figs.pfc-d)] at strong anisotropy 
ft = 0.2. The effect of j c ^ Q is considerable. Even at 
7 C = 0.2 < (7||), relaxation of the central spin is already 
much faster than at 7 C = 0. For 7 C > 7|| , we find qualita- 
tive changes: the local minimum disappears, relaxation 
becomes almost complete and furthermore becomes ex- 
ponential [see Fig. [2^b)] . At 7 C > 8, the exponential 
relaxation rate saturates at a value that does not depend 
on 7c any longer. Figures [2^c-d) show analogous results 
for the nuclear spin bath with the Hamiltoninan ([9| and 
1 = 1 with 7g = 47*, which corresponds to the same 
characteristic splitting of energy levels by QC. It shows 
that by changing the size of spins and form of the cou- 
pling but keeping the same characteristic 7 C , the form of 
the central spin correlator does not change. 

To better understand the change of behavior with 
growing j c , it is instructive to look at the dynamics of 




FIG. 3: For f3 = 0.2: (a) Typical Overhauser field dynamics 
for 7c = (pink line) and 7 C = 8 (blue), (b) Real time 
Overhauser field correlator C^ tll {t) = (B nz (t)B nz (0)). 



B nzi the Overhauser field component along the z-axis. 
Figure [3] shows examples of B nz (t) starting from a ran- 
dom initial condition for all spins. When 7 C = 0, B nz 
is practically frozen. However, for 7 C > 7|| values, B nz 
quickly fluctuates with the amplitude of the typical Over- 
hauser field strength. Figure [3jb) shows that in the latter 
case, the bath spin correlator C2 ath (£) = (B nz (t)B nz (0)) 
decays during t < l/^f c to a smaller but nonzero value. 
For 7c ^> 71 1 , nuclear spins simply precess around their 
local static fields from the QC. Fluctuations that are 
seen in Fig. |3ja) are then merely due to the difference 
of precession frequencies and precession axis directions 
for different nuclear spins in one dot. In contrast, when 
7c/7|| < 1, nuclear spin precessions are synchronized by a 
stronger hyperfine coupling that suppresses fluctuations 
of B nz . 

When 7c > 7||, there can be two distinct regimes of 
central spin polarization dynamics. The first regime ap- 
pears when fluctuations of the Overhauser field are so fast 
that the adiabaticity conditions break down and the cen- 
tral spin polarization cannot follow the direction of the 
Overhauser field. This most likely can happen when B nz 
passes through zero values and the Landau- Zener tran- 
sition probability, plz — 1 — exp(— 7r(B n ±) 2 /v) is sub- 
stantially different from unity. Here v = (dB nz /dt)B nz =o 
and B n ± ~ j_\_y/N is the typical value of the Overhauser 
field transverse to the z-axis direction. In this case, each 
time B nz changes sign, the central spin has substantial 
probability of not following the Overhauser field so that 
its dynamics become stochastic with exponential relax- 
ation of the central spin correlator 23 . To estimate plz 
we note that, according to Fig. |3^b), when 7 C > 7||, the 
time l/7c sets the scale for the correlator decay time of 
the Overhauser field. The latter changes during this time 
by the amount SB nz ~ ^y\\y/N. Hence, the rate of change 
of the Overhauser field is v ~ j\\y/Nj c ^ and exponential 
relaxation can occur when two conditions are satisfied: 



(3) 



7c>7||, and rj = ft 2 ^ Vn/j c < 1. 

For our numerical test with TV = 700 and f3 = 0.2, we find 
that (19) is satisfied when J c /j\\ ^ 1- This result is in 
agreement with Fig. [2] The exponential relaxation time, 
r re i, roughly corresponds to the value of l/j c a t which 
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FIG. 4: (a) Real time spin correlator and (b) frequency power 
spectrum P(u) = f dte tLOt (a z (t)a z (0)) for different values of 
external out-of-plane magnetic field B z ; j c = 8, and f3 = 0.2. 




FIG. 5: Central spin correlator at f3 — 1 in zero external field: 
(a) Real time correlator and (b) frequency power spectrum. 



Plz ~ 1/2, i.e. r re i ~ l/[/3 2 7|| V^V]. For the hole-doped 
dot pl22| we assume TV = 10 5 , (3 = 0.2, 7c /7|| = 25, 
which gives 77 = 0.4 < 1, i.e. it agrees with the observed 
Lorentzian shape of the hole spin noise power spectrum 
in [10]. Considering that corresponds to several 

microseconds in hole-doped dots, we find the relaxation 
time to be a fraction of a microsecond, which also agrees 
with the experimentally measured value r re \ ~ 0.4 [is at 
a zero external fields Our model is also in good agree- 
ment with other experimental observations: For example, 
when an external out-of-plane magnetic field was applied, 
the central spin relaxation was suppresse d 10 * 11 * when this 
field exceeded j\\y/N. In Fig. J^a) we confirm this fact 
numerically. Figure [4jb) also shows our numerical results 
for the effect of an applied magnetic field on the hole 
spin noise power spectrum, P(u) = j dte zu;t (a z (t)d- z (0)), 



which is in good agreement with experimental measure- 
ments of this spectrum in external fieldf^^R 

The second regime corresponds to the case when fluc- 
tuations of the Overhauser field are strong but the central 
spin follows the direction of the Overhauser field adiabat- 
ically. This happens when 

7c > 7|| » 7 1 >1 - ( 4 ) 

For electron- doped InGaAs dots, condition Q would 
likely be satisfied because of a lack of anisotropy (/? = 1). 
Assuming that such a dot has N ~ 10 5 nuclei and 
7c/7|| = 3, we find n ~ 10 2 ^> 1, i.e. the central spin 
dynamics is well within the adiabatic regime. For such 
conditions, the central correlator has to follow the corre- 
lation pattern of the Overhauser field, as in Fig. j3^b). 

Figure [5|a) shows our results for the real time corre- 
lator for j3 = 1 (electron-doped dots) and N = 700. The 
case 7 C = 3 corresponds to conditions Q . The first min- 
imum of C2(t) in Fig. [5|a) is due to dephasing effects. 
Note that it is not destroyed by Overhauser field fluctu- 
ations, unlike the case with j3 = 0.2. At longer times, 
C2(t) qualitatively follows the Overhauser field correla- 
tion pattern, i. e. it decays to a small but non-zero value 
during a time ~ l/7c? followed by a long relaxation tail. 
Figure [5|b) shows that a specific feature of the regime 
Q , which distinguishes it from the case with j c < 7| | , is 
trie appearance of a shoulder in the low frequency peak of 
the spin noise. An additional feature of the power spec- 
trum at n > 1 is the presence of a second broad small 
amplitude Gaussian peak at high frequencies. 

In conclusion, we identified three regimes with distinct 
central spin dynamics in the presence of QC at low tem- 
peratures and weak external fields: (i) the regime of ex- 
ponential relaxation of the spin correlator, which is de- 
fined by Eq. ( [19] ) ; (ii) the regime with the central spin 
following Overhauser field adiabatically, which is defined 
by Eq. Q; and (hi) the regime of weak QC, < 7 C < 7||, 
which is qualitatively similar to j c = 0. We showed that 
hole-doped InGaAs dots^l likely correspond to the expo- 
nential relaxation regime and that electron-doped dots 
correspond to the regime (ii). Regime (iii) is potentially 
applicable to electrostatically defined dots with a nearly 
perfect atomic lattice. 
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Supplementary Material 

In this Supplementary we discuss our numerical ap- 
proaches in more detail and provide representative results 
for numerical simulations of the central spin problem for 
larger spin bath sizes and higher spins (with / = 1). 
We discuss our alternative numerical approach, which is 
based on solving classical Landau-Lifshitz equations with 
quadratic anisotropics for nuclear spins. We show that 
although physically justified changes of the model Hamil- 
tonian may lead to observable quantitative differences, 
our basic conclusions about the role and effects of the 
quadrupole couplings are robust. 



I. NUMERICAL FRAMEWORKS FOR 
MODELING CENTRAL SPIN RELAXATION 

A. TDMF framework to simulate static field model 
of a spin-1/2 spin bath 

In this model, we assume that all spins have size 1/2 
and then we mimic the effect of QC by introducing ran- 
dom quenched magnetic fields acting on each nuclear 
spin, so that the full effective Hamiltonian for numeri- 
cal simulations reads: 



N 

H = Ba° + + 

+ 7±(*2*x + W +'£(**«')). (5) 

where the index i = stands for the central spin, a l a is 
the a-component of the Pauli operator for the i-th nu- 
clear spin. The vector rt % points in a random direction 
and 7* is the strength of the quenched field acting on 
the z-th nuclear spin. We will assume that 7* have the 
uniform distribution in the interval (0,27 c ), i.e. with 
mean (7*) = j c . Formally, the substitution of the QC 
by the random quenched magnetic field changes sym- 
metry properties of this coupling. However, we believe 
that the Hamiltonian ([5| captures the basic qualitative 
semiclassical picture that leads to the enhanced central 
spin relaxation. The main advantage of this approach 
is that there is an efficient well-scalable numerical algo- 
rithm, based on the time-dependent mean field (TDMF) 
dynamics, to simulate the quantum mechanical evolution 
of the central spin relaxation^. Although it is not exact 
by construction, the TDMF algorithm showed results for 
the central spin relaxation that are indistinguishable from 
exact quantum mechanical simulations when the number 
of nuclear spins is relatively large N > 20 and when the 
evolution is averaged over a random nuclear spin distri- 
bution at high temperature of the nuclear spin batrW An 
additional advantage of this model is that the evolution 
of spins- 1/2 in a magnetic field can be relatively quickly 
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computed so that it is possible to simulate the dynamics 
with a sufficiently large size of the nuclear spin bath. 

Following the TDMF theory, we approximate the state 
vector of the total system as a product = 
\ u o) T\i=i \ u k) °f the single-spin vectors \uj), (j = 
0, . . . N). Then, at each time step we update the state of 
each spin by considering its evolution with the effective 
Hamiltonian 



H i eR = h i {t)a i , i = 0,l,...,N 



(6) 



where the effective field hi(t) acting on the i-th spin is 
calculated, at each step, according to 



hn = B 



N 
i=l 



7||^* + 7xK* + ^v), (7) 



hi = 7K* + 7i(0£A + *i»)+7; 



1>\ 



(8) 



where we defined cr 3 (t) = Tr[p(£)cr J ], and where p(t) is 
the density matrix that corresponds to the pure state 

*(*). 

More information about the theoretical justification of 
this approach can be found in Ref. pQ. Although this 
reference does not discuss effects of additional quenched 
fields, one can easily trace that the same steps lead to 
when all spins have a value 1/2. 



B. TDMF framework to simulate spin bath with 
spins I > 1/2 with quadrupole coupling 

To simulate dynamics with the Hamiltonian 

N 

H = Y1 (itftzSz + 7±(4A + liySy)) + 9zB z S z + 
i=l 



(9) 



+9: 



N i 

,B X S X + g v B y S v + J2~Y ' ni)2 ) ' 



with / > 1/2, we make a similar mean field approxima- 
tion that the central spin moves in the effective time- 
dependent Overhauser field: 



H° s = h (t)&°, 



(10) 



where 



N 



ho = B + ^ l r z z + f ± (i l x x + i l y y)- (ii) 

i=l 

There is, however, a complication for nuclear spins be- 
cause in the spirit of the TDMF- approach, only the cou- 
pling to the central spin should be treated in the mean 
field approximation, while all local fields should be in- 
cluded exactly in numerical simulations. Therefore, the 
effective Hamiltonian for the nuclear spins is given by 



Hiv = H i c {t) + H i Q = K{t)-r 



((Ii ■ n,) 2 ) , 



7Q 

2 eff = -"cW t"q = "»W ' 1 "T- -y 

where i = 1,2, . . . , N, and where 

h i = 1 \ f ° z z + 1 i L {alx + a° y y). (12) 

We split the effective Hamiltonian for the z-th nuclear 
spin into the part H % c that is responsible for interaction 
with the central spin and the part Hq that is responsible 
for the quadrupole coupling. The reason for this splitting 
is that, in numerical simulations, we should update the 
state vector by the evolution operator 



U 



iHl ff dt 



(13) 



at each time step dt. For spins I > 1/2 this matrix 
exponent cannot be written explicitly in a compact form. 
Instead, to perform simulations, e^g. with / = 1, we used 
the Suzuki- Trotter decomposition^ 



U 



e iH % Q dt/2 e iHldt e iH % Q dt/2 



(14) 



which is the unitary approximation that is correct up to 
the order 0(dt 2 ). For 1 = 1 and I = 3/2, all matrix 
exponents in ( 14 ) can be readily written explicitly. 



C. Spin bath model with classical Landau-Lifshitz 
dynamics 

For spins I > 3/2, matrix exponent of Hq is not 
available in analytical form, however, one can apply an- 
other approximation U « (l+z^ ff dt/2)(l-^ ff dt/2)- 1 
which is again unitary and exact up to O(dt). It reduces 
the calculation of the evolution operator to calculations 
of the matrix inversion which is readily available. In- 
stead of this approach, however, we explored the case of 
spin bath dynamics with large nuclear spins by assuming 
a classical approximation of Landau-Lifshitz dynamics 
with the Hamiltonian ([9|. 

We considered the classical spin Hamiltonian, which 
correctly takes into account the form of the quadrupole 
coupling: 



H = BS 



N 

£ 

i=i 



7q , 



i\2] 



(15) 



leading to a standard dissipationless Landau-Lifshitz dy- 
namics: 



hj x Sj, S = Hq x S, 



(16) 



N 



h = B + 5^7f|4« + 7 i(4* + 4y), (17) 
hi = 5 f | ^z + 7 i(^ + 5 y y)+7^(s i n i K. (18) 
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FIG. 6: The real time spin-spin correlator for nuclear spins 
baths with the Hamiltonian |5| and N = 10000 nuclear spins; 
f3 = 0.2, and B = 0. Different curves correspond to different 
values of j c . Time is in units of Averaging is over 5000 

runs with different random initial conditions. 



The classical treatment of nuclear spins should be ad- 
ditionally justified in our case by large sizes of nuclear 
spins (9/2 for indium isotopes). For simplicity, we will 
normalize to one all spin sizes in this classical model. 
The advantage of the classical approach is its speed. It 
also correctly preserves all symmetries of the couplings. 
Its disadvantage is its inability to take rigorously into ac- 
count the discreteness of spin states, which turned out to 
be not an important factor to obtain qualitative behavior 
of the central spin problem. 



II. SCALING OF CENTRAL SPIN 
CORRELATORS WITH THE NUMBER N OF 
NUCLEAR SPINS 

As an additional test of our theory, we performed a 
number of simulations of the central spin dynamics with 
different sizes of the nuclear spin bath in the range from 
N = 250 to N = 10000 nuclear spins. To obtain results 
with bigger N we enhanced precision by decreasing the 
size of the time step in simulations to dt = 0.00005 and 
reducing the statistics of averaging over different random 
initial conditions from 20000 runs for N < 1000 to 5000 
runs at N = 10000. In the main text, we made two 
analytical predictions which are particularly suitable to 
test by changing N. First, we predict that the transition 
between the exponential regime and the regime at which 
the central spin polarization follows the Overhauser field 
adiabatically happens in the region of parameter values, 
at which the following condition is satisfied: 



with the exponential relaxation regime corresponding to 
r] < 1. 

Figure [6]^a) shows the effect of increasing 7 C for spin 
bath with N = 10000 nuclear spins. At 7|| =1 and j3 — 
0.2, condition ( p~9| ) would be satisfied for j c « 4. Indeed, 
Fig. |6|a) shows that near the value of j c = 3 the local 
dephasing minimum at short time-scales is still present 
but almost disappeared. At j c = 12 the relaxation is 
exponential. One can compare this figure with analogous 
Fig. 2 in the main text for N — 700. In the latter case, 
similar disappearance of the dephasing minimum of C2 (t) 
was observed at a smaller value of 7^ ~ 1, which is in 
agreement with the \fN factor in (19). 



Another apparent difference of cases with N = 700 and 
N = 10000 is that at the larger N the local dephasing 
minimum of C2(t) shifts toward the shorter time scales. 
This is expected because this minimum corresponds to 
the typical rotation time of the central spin around the 
Overhauser field. This field increases with N as y/~N, 
which explains the shift of the local minimum. We note 
also that increasing N does not change our conclusion 
about the absence of substantial relaxation of the central 
spin at 7 C = during times longer than 1/7||. 

Our second analytical prediction states that in the 
exponential regime of relaxation, the relaxation rate is 
given by 



Trel 



(20) 



where (3 is the anisotropy factor, 7|| is the typical strength 
of hyperfine coupling and N is the number of spins. Co- 
efficient A is a numerical factor, which is not universal 
and may depend on the specific characteristics of the nu- 
clear spins, for example the sizes of the nuclear spins and 
the distribution of parameters. Nevertheless, considering 
the generality of the arguments leading to Eq. (20), this 
coefficient is expected to be of the order unity. 



In order to test prediction (20 ), we considered a nuclear 
spins- 1/2 bath with anisotropy /3 = 0.2 and 7 C /7|| = 12. 
According to (19) the condition for exponential relax- 

12 and /3 0.2 for 
90000 nuclear spins. 



ation, 77 < 1, is satisfied at j c - 
nuclear spin baths with up to N 

Figure [7] shows our numerical results for the relaxation 
times of correlators that were obtained for different TV 
by the best exponential fits of the correlator relaxation 
curves (blue boxes). Red curve shows the best fit of these 
numerically obtained points by Eq. (20) with a single 



77 = 2 7||ViV/7c~l, 



(19) 



free fitting parameter A. The best fit was obtained for 
A « 0.45, which is of the order unity, as expected. Fig- 
ure [7] demonstrates quite a reasonable agreement with 
prediction (20). 
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FIG. 7: Relaxation time as a function of the number N of nu- 
clear spins-1/2 in the bath. 7|| = 1, ft = 0.2, j c = 12, which 
corresponds to the exponential regime. Time is in units of 
1/7||. Blue boxes correspond to results of nu mer ical simu- 
lations and dashed red curve is given by Eq. (20) with the 
numerical factor A 



■ 0.45. 



III. CENTRAL SPIN CORRELATORS WITHIN 
TDMF APPROACH WITH SPIN-1 NUCLEAR 
SPINS 

A. Relaxation for different j c — 7q/4 

To show that our results are not specific for nuclear 
spins-1/2, we performed TDMF simulations for spins- 1, 
treating the quadrupole coupling exactly with 2nd order 
Suzuki- Trotter decomposition^, as explained in Section 
LB. 

In order to make quantitative comparison with results 
for spin- 1/2, we note that the value 2j c was the split- 
ting of energy levels in the spin- 1/2 Hamiltonian, which 
is induced by static fields. The corresponding term in 
the Hamiltonian with quadrupole coupling, (jq /2){n-I) 
introduces the level splitting of the size 7q/2 between 
nearest energy levels. Hence, for \I\ = 1 we introduce 
parameter j c = 7q/4, which characterizes the level split- 
ting by quadrupole field. Since we expect that j c is the 
main parameter that determines the effect of quadrupole 
coupling, the behavior of spin correlators at same values 
of 7 C should be similar for spin- 1/2 and spin-1 nuclear 
spin bathes. 

Figure [8] shows the results of our simulations for spin- 
1 nuclear spin bath. It reproduces Figs. 2(c,d) of the 
main text. Distributions of parameters Y\\ an d J± were 
chosen the same as in the case of spin- 1/2 bath in the 
main text. Quadrupole coupling for i-th nuclear spin 
was defined 7q = 47*, where 7* was chosen according to 
the same distribution as in the main text for the spin- 1/2 
bath. 
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FIG. 8: (a) The real time spin-spin correlator (blue curve) 
for nuclear spin I — 1 baths with different magnitudes of 
the quadrupolar coupling. We introduce a characteristic level 
splitting parameter j c = /4 in oder to compare predictions 
of spin-1/2 and spin-1 baths; f3 = 0.2, N = 700, and B = 0. 
(b) Exponential fit (dashed red) of the spin correlator at j c = 
12 in units of 7|| (red). Time is in units of 1/7||. 



Comparing Figs. 2(c,d) and Figs. 2(a,b) of the main 
text, which describe, respectively spin-1 and spin-1/2 
baths, we find that the differences between them are 
not essential despite difference of the nuclear spin sizes 
and the form of the considered nuclear spin coupling to 
strain fields. As in the analogous example of spin-1/2, 
the transition to the exponential regime happens at a 
characteristic parameter value j c ~ 1, with relaxation 
exponent saturating near j c ~ 8. More detailed exami- 
nation reveals small differences, for example, for 7 C = 12, 
the characteristic relaxation time obtained by the best 
exponential fit for spin-1 bath is equal to t^j 1 ~ 0.50, 
while for the spin-1/2 bath in Fig. 2 of the main text it 

was t^ 1 ^ 2 ~ 0.39. Such quantitative differences up to a 
coefficient of order unity are expected and do not affect 
any of the basic predictions of our theory. 



B. Relaxation at fixed 7q but different N 

Figure [9] tests the 1 / prediction for the relaxation 
time in the exponential regime for spin baths with nu- 
clear spins 1 = 1 and N in the range from N = 200 
to N = 4000. Again, as in the case of spins-1/2, the 
agreement with theoretical prediction is quite reasonable. 
Quantitative differences between cases with defferent nu- 
clear spin sizes mainly reduce to a small difference of the 
corresponding coefficients A: A w 0.56 for spins-1 bath 
and A w 0.45 for spins-1/2 bath. 



C. Effect of coupling randomness 

Here we use the TDMF approach for / = 1 to inves- 
tigate the role of coupling randomness. In all our previ- 
ous examples, we assumed that parameters 7^, j]_ and 

7* have a broad distribution with the variance close to 
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FIG. 9: Relaxation time as function of the number N of 
nuclear spins with I = 1 in the bath. ^ = 1, /3 = 0.2, 
7c = 7q/4 = 12, which corresponds to the exponential regime. 
Time is in units of 1/7||. Blue boxes correspond to results 
of numerical simulations and dashed red curve is given by 
Eq. (20) with the numerical factor A « 0.56. Averaging is 



performed over 5000 runs with different initial random con- 
ditions for all spins. 




FIG. 10: Spin-spin correlators in the absence of randomness 
for parameters 7^ , j]_ and 7*. The nuclear spin bath has I — 1 
with random choice of anisotropy vectors n\ Level splitting 
parameter is -y c 7q/4; /3 0.2, N 700, and B = 0. 
Averaging is over 10000 runs with random initial conditions. 



the mean value. However, most of our discussion for the 
transition to exponential regime and the relaxation time 
value did not involve this randomness assumption. 

We performed simulations, analogous to the ones in 
Fig. |8j but assuming that 7^ = 7| | , 7j_ = 7± and 7* = j c 
for any i = 1, ... TV. We kept the assumption of random 
n l only. While such a spin bath may be considered un- 
physical, it is instructive to look at its effects in order to 
illustrate the effect of parameter randomness, 



7 C ~ 1 remains true, and at 7 C = 8 we find an almost ex- 
ponential relaxation, as usual. If we compare relaxation 
rates, we find that they are not influenced by parameter 
randomness. For example, for the case in Fig. [8] the best 
exponential fit of j c = 8 curve corresponds to r re \ ~ 0.50 
and for the case in Fig. [To] we found r re i ~ 0.51. This 
example illustrates the universality of behavior in the 
exponential regime if it is achieved. While for small and 
intermediate j c values, strong changes of parameter dis- 
tribution, in principle lead to noticeable changes of re- 
laxation pattern, after reaching the exponential regime 
the differences between cases with different distributions 
of parameters become practically unnoticeable. 

At the regime with 7 C /7|| < 0.5, the above behavior 
also changes insignificantly. The difference becomes clear 
if we look at the 7 C = case. Then the model with 
identical parameters predicts no relaxation at all after 
going through the dephasing minimum. 

Considerable changes of the behavior are observed at 
7c/7|| °° 1 5 for which a very complex relaxation pattern is 



shown in Fig. 10 This behavior is not surprising either. 
Indeed, this case corresponds to a situation at which fluc- 
tuations of the Overhauser field are substantial but the 
central spin follows those fluctuations almost adiabati- 
cally. For constant parameters in this regime, the corre- 
lator of the Overhauser field can show unusually large os- 
cillations, which at normal conditions would be smeared 
by parameter randomness. The central spin follows that 
Overhauser field and, in this regime, its correlator re- 
flects oscillating behavior of the Overhauser field. This 
explains the unusual relaxation at j c /j\\ ~ 1 m Fig. 
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IV. CENTRAL SPIN CORRELATORS IN 
SEMICLASSICAL LANDAU-LIFSHITZ 
DYNAMICS 

The complexity of simulating spin baths with higher 
sizes of nuclear spins increases. To get an insight in 
spin bath effects with large values of spins, we simu- 
lated the nuclear spin bath dynamics classically within 
the Hamiltonian (15). As it was expected from the semi- 



classical character of our theory, we obtained essentially 
the same relaxation behavior of the central spin as in the 
case of the quantum mechanical TDMF framework with 
quenched magnetic fields. The quantitative differences 
should be attributed to normalization [C^O) = 1/3 for 
classical unit spin rather than 6*2(0) = 1 for quantum 
spin- 1/2] and to different functional forms of the static 
field terms assumed in the two frameworks. Since all re- 
sults that we obtained are qualitatively the same as in 
the case of the TDMF framework with quenched mag- 
netic fields, we will only show a few examples. 



Figure 10 shows results of such simulations. First, we 
note that the transition to the exponential regime above 



We performed simulations of the model (15) for a 
central spin coupled to N = 700, and in a few in- 
stances to N = 4000 nuclear spins. We chose quadrupole 
couplings to be uniformly distributed random numbers 
7q = 7Q * RNDi, where [RNDi e (0, 2)], and the direc- 
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FIG. 11: Real time spin-spin correlator produced by classical 
Landau-Lifshitz dynamics with Hamiltonian ( |15| ). Here N = 
700, (7| | ) = 1, P = 0.3. Results are averaged over 20000 
runs with random initial conditions for central and nuclear 
spins. The correlator of classical spin with \S\ = 1 starts at 

e 2 (o) = i/3. 
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FIG. 12: Real time spin correlator in classical Landau-Lifshitz 
dynamics with N = 4000 nuclear spins, (71 1) = 1, /3 = 0.2. 
Results are averaged over 5000 runs with random initial con- 
ditions for central and nuclear spins. 
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FIG. 13: Real time spin-spin correlator C2(t) = (S z (t)S z (0)) 
in c lassical Landau-Lifshitz dynamics with the Hamiltonian 
|l5| ), at 7q = 0, up to times (a) t = 5/ (711) and (b) t = 
0.5/ (7| 1 ) • Different curves correspond to different sizes of the 
coupling anisotropy, f3 = (7_l)/(7||). Here B = 0, N = 4000, 
(7| I ) = 1. Averaging is (a) over 5000 and (b) over 25000 
simulations with different random initial states of central and 
nuclear spins, as well as hyperfme couplings. 



EFFECT OF NUCLEAR STATIC FIELD 
ANISOTROPY 



tion of the anisotropy axis rti was chosen randomly for 
each nuclear spin. 

According to Fig. [TTJ as in the case of quenched fields, 
there is a critical value of 7q above which the central 
spin relaxation becomes almost exponential and the lo- 
cal minimum due to the coherent central spin precession 
disappears. To show that this behavior is not specific for 
N = 700 nuclear spins, Fig. 12 shows C2(t) for a larger 



number of nuclear spins (N = 4000). 



Finally, Fig. 13 shows the numerically calculated spin 
correlator for the model with the Hamiltonian (15) at 
7q = with N = 4000. Different curves correspond 
to different values of the coupling anisotropy j3. Fig- 
ure 13 [b) resolves the part of Fig. 13 'a) with t < 0.5. 



In the main text, we assumed that quenched fields 
can equally probably point in any direction. In real- 
ity, we expect that anisotropy of the lattice is reflected 
in anisotropy of quadrupole coupling fields. To explore 
this case within TDMF approach, we considered a cubic 
anisotropy by choosing random quenched fields 7 C * in 
spin- 1/2 bath having independent random components 
7cx' lcy> an d llzi eacn uniformly distributed in the in- 
terval ("7c, 7 C ). Such random field vectors are uniformly 
distributed inside a cube rather than a sphere. Figure [14} 
which is the counterpart of Fig. 3 in the main text, 
shows that the correlation function of the Overhauser 
field [Fig.[l4jb)] is quite different from the one in Fig. 3 of 
the main text but the sample of a fluctuating trajectory, 
shown in Fig. [l4[a) is similar to the isotropic case. This 
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FIG. 14: For f3 = 0.2 and cubic anisotropy of 7 C : (a) Typical 
Overhauser field dynamics for j c = (pink) and j c = 8 (blue) ; 
(b) real time Overhauser field correlator. N — 700. 



means that our basic discussion about the breakdown 
of adiabaticity should equally apply to the case of cubic 
anisotropy, leading to the same quantitative predictions. 
Indeed, Fig. 15 , which is the analog of Fig. 2(a,b) in the 
main text, shows that a cubic anisotropy of quenched 
fields does not affect our conclusions about relaxation 
rates at different strengths of j c . Effects of additional 
in-plane or out-of-plane anisotropy of j c can be more 
pronounced because in-plane anisotropy favors stronger 
fluctuations of B nzi while out-of-plane anisotropy sup- 
presses fluctuations of B nz . 
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FIG. 15: Spin correlators at cubic anisotropy of 7 C : (a) Real 
time spin correlator for different magnitudes of the quenched 
field 7c, f3 = 0.2, B = 0; (b) exponential fit (dashed black) of 
spin correlator for j c = 8 (blue). N — 700. 
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